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ABSTRACT 

Massive early-type galaxies have undergone dramatic structural evolution over the last 10 Gyr A compan- 
ion paper shows that nearby elliptical galaxies with M» > 1.3 x IO^'M© generically contain three photometric 
subcomponents: a compact inner component with effective radius Re <^i kpc, an intermediate-scale middle com- 
ponent with Re ~ 2.5 kpc, and an extended outer envelope with Rg « 10 kpc. Here we attempt to relate these 
substructures with the properties of early-type galaxies observed at higher redshifts. We find that a hypothetical 
structure formed from combining the inner plus the middle components of local ellipticals follows a strikingly 
tight stellar mass-size relation, one that resembles the distribution of early-type galaxies at z « 1. Outside of the 
central kpc, the median stellar mass surface density profiles of this composite structure agree closest with those 
of massive galaxies that have similar cumulative number density at 1.5 < z < 2.0 within the uncertainty. We 
propose that the central substructures in nearby ellipticals are the evolutionary descendants of the "red nuggets" 
formed under highly dissipative ("wet") conditions at high redshifts, as envisioned in the initial stages of the 
two-phase formation scenario recently advocated for massive galaxies. Subsequent accretion, plausibly through 
dissipationless ("dry") minor mergers, builds the outer regions of the galaxy identified as the outer envelope in 
our decomposition. The large scatter exhibited by this component on the stellar mass-size plane testifies to the 
stochastic nature of the accretion events. 
Subject headings: galaxies: photometry — galaxies: structure — galaxies: evolution 



1. INTRODUCTION 

Recent observations have established that high-redshift 
early-type galaxies (ETGs) are more compact (Daddi et al. 
2005; Trujillo et al. 2006; Damjanov et al. 2011) and have 
higher stellar velocity dispersions (Cappellari et al. 2009; On- 
odera et al. 2009) than nearby ETGs of the same stellar mass. 
Since z « 2.5, "red nuggets" on average have doubled in stel- 
lar mass and increased their size by a factor of 3 - 4 (Buitrago 
et al. 2008; van Dokkum et al. 2010; Papovich et al. 2012; 
Szomoru et al. 2012). Some appear to have a disk-like mor- 
phology (van der Wei et al. 201 1), and they are always bluer in 
the outskirts (Gargiulo et al. 2012). The accumulated evidence 
suggests that massive ETGs build up inside-out through non- 
dissipational processes (Bezanson et al. 2009; van Dokkum et 
al. 2010). 

These developments severely challenge classical models of 
elliptical (E) galaxy formation such as monolithic collapse 
(Larson 1975) and binary major mergers (Toomre & Toomre 
1972; Negroponte & White 1983). Instead, the current evidence 
favors a "two-phase" scenario (Oser et al. 2010; Johansson et 
al. 2012). Intense dissipational processes such as cold accretion 
(Dekel et al. 2009) or gas-rich mergers rapidly build up an ini- 
tially compact progenitor After star formation is quenched, a 
second phase of slower, more protracted evolution is dominated 
by non-dissipational processes such as dry, minor mergers. 

For such a scenario to work, the rate of minor mergers (Bluck 
et al. 2012; Newman et al. 2012) has to be consistent with the 



prevalence of faint companion galaxies at high redshifts. Per- 
haps other mechanisms may also be important, such as major 
merger (Bemardi et al. 2011; Prieto et al. 2013) or AGN- 
induced expansion (Fan et al. 2008). Given that large uncer- 
tainties remain in identifying the relevant physics and quantify- 
ing their detailed balance, a key sanity check is to see whether 
one can separate out the fossil bodies of compact z > 1 .5 galax- 
ies at the core of nearby galaxies. This separation is success- 
ful only when several clear and solid predictions, especially the 
similarity between the core components and the high-z compact 
progenitors in mass versus size relation and detailed mass pro- 
file, are met. These stringent predictions are the basis of this 
study to verify or falsify. We assume a ACDM cosmology with 
n^ = 0.27, flA = 0.73, and Hq = 73 kms"' Mpc-\ 

2. OBSERVATIONAL MATERIAL 

This work uses the three-component models of nearby Es 
from the Carnegie-Irvine Galaxy Survey (CGS; Ho et al. 201 1), 
a systematic study of 605 bright {Br < 12.9 mag) southern 
(S < 0°) galaxies. The currently completed photometric part 
of CGS includes BVRI images obtained using the 100-inch du 
Pont telescope at Las Campanas Observatory, using a 8f9x8f9 
detector with a pixel scale of 0."259, under ■^l" seeing condi- 
tions. More details on the observations and data reductions can 
be found in Ho et al. (201 1) and Li et al. (201 1). 



Huang et al. (2013) apply GALF I T (Peng et al. 2010) to per- 
form a detailed two-dimensional analysis of the V-band images 
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Fig. 1. — The redshift evolution of size and its relation with stellar mass for nearby Es and high-z ETGs. The left side of panel (a) summaiizes the size 
distributions of single-component models (black circles) and the inner (red triangles), middle (green squares), and outer (blue pentagons) components of the three- 
component models of CGS Es. The median of each distribution is marked by a horizontal bar. The right side illustrates the redshift evolution for high-z ETGs 
with M» > IO'^'Mq (grey open circles). The solid and dashed lines give the running median and its associated first and third quartile through five redshift bins 
(0.50 < z < 0.75; 0.75 < z < 1.0; 1.0 < z < 1.4; 1.4 < z < 2.0; 2.0 < z < 2.7). Panel (b) shows the mass-size relation for high-z ETGs, separated into those with 
z < 1.5 (grey stars) and z > 15 (grey skeletal triangles). The 21 massive ETGs at 1.5 < z < 2.5 from Szomoru et al. (2012) are marked with cyan circles. Each of 
the three subcomponents of the CGS sample is plotted separately, as well as for the inner-Fmiddle components combined (orange diamonds). Also highlighted are 
the median properties for the inner-Fmiddle components of the CGS Es and for the massive galaxies at 1 .0 < z < 1 5 and 1 .5 < z < 2.0 selected by their cumulative 
number density from Patel et al. (2013). Typical eiTor bars for the subcomponents are shown on the lower-right comer. The best-fit relation of Guo et al. (2009) for 
local SDSS ETGs is plotted with a dashed line. For consistency with the high-z measurements, the effective radii used here are the so-called mean, or circulaiized, 
radii, such that Re,c = R^ \/b/a, where b/a is the axis ratio. 



of 94 CGS Es. Their decomposition reveals that, contrary to 
popular perception, the global light distribution of most nearby 
Es have more complicated structures than can be described by 
a single Sersic (1968) profile. Instead, 70 out of the 94 Es can 
be fit with three distinct subcomponents: (1) an inner, compact 
center with effective radius Re <^^ kpc comprising a small frac- 
tion of the light (/ « 0.1-0.15); (2) an intermediate-scale, mid- 
dle structure with R^ « 2.5 kpc and / w 0.2-0.25; and (3) an 
outer, dominant, extended, moderately flattened envelope with 
7^, w 10 kpc and/ « 0.6. 

The image decomposition was performed without ascribing 
any physical significance to the substructure. Nevertheless, 
the fact that the individual subcomponents define distinct se- 
quences on the size-luminosity plane, and that three-component 
fits often seemed naturally preferred, hint at the possibility that 
the subcomponents reflect underlying physical reality. Here we 
compare some key properties of the three-component models 
of local Es with observations of high-z ETGs to shed new light 
on the structural evolution of massive galaxies. 

For this purpose, we compiled stellar masses and effective 
radii for 1323 high-z ETGs with stellar masses M, > IO'^^Mq, 
as follows: 352 (0.2 < z < 2.7) from Damjanov et al. (2011); 
910 (0.2 < z < 1.2) from COSMOS (ScoviUe et al. 2007; us- 
ing morphological classification from Scarlata et al. 2007 and 
Tasca et al. 2009); 32 (1.5 < z < 3, choosing Sersic indices 
n > 2) from the GOODS-NICMOS survey (Conselice et al. 
2011); 8 (z w 1.6, n > 2) from Papovich et al. (2012); and 21 
(1.5 < z < 2.5) from CANCELS (Grogin et al. 2011, Koeke- 
moer et al. 2011), as analyzed by Szomoru et al. (2012). Al- 
though the sample is in many respects heterogeneous (in terms 
of observed bandpass, selection criteria, modeling method), the 
general trends we explore in this paper should not be strongly 



affected, as demonstrated by Damjanov et al. (201 1). 

3. LINKING NEARBY AND HIGH-Z EARLY-TYPE GALAXIES 

Figure 1(a) illustrates the redshift evolution of galaxy size 
for massive ETGs. Although the scatter is substantial, the size 
increase with decreasing z is confirmed (e.g., Damjanov et al. 
2011; Cimatti et al. 2012). The median circularized effective 
radius increases from 0.8 kpc at 2.0 < z < 2.7, to Rg.c ~ 1 .9 kpc 
at 1 .4 < z < 2.0, to /;,,, « 3 kpc at 0.50 < z < 0.75. By z = 0, Es 
have reached a median Rec = 7.5 kpc according to the single- 
component fits applied to the subsample of 35 massive CGS Es 
that satisfy M, > 10" M© (soHd black points). The CGS mass 
cut roughly corresponds to twice the mass limit of the high-z 
sample. 

How are the subcomponents of nearby Es related to the high- 
z systems? Judging by the distribution of sizes for the subcom- 
ponents shown on the left side of Figure 1(a), the typical effec- 
tive radius of the inner component (median Re^ = 0.9 kpc) does 
overlap substantially with that of high-z red nuggets; however, 
their stellar masses, only ^10%-15% of the total, are about a 
factor of 3-5 lower than the high-z systems (Figure lb). The 
middle component, with median R^.c = 3.7 kpc, matches well 
the sizes and masses of ETGs at z ~ 0.5, but they are still on 
average more extended than those at z > 1 .0. On the other hand, 
with median R^.c =17.5 kpc, the outer component of local Es is 
clearly larger than the globally averaged size of most ETGs at 
any redshift. 

The mass-size relation offers another vantage point. As de- 
scribed in Huang et al. (2013), single-component fits to the 
CGS Es produce sizes that follow well the mass-size relation of 
local ETGs selected from the Sloan Digital Sky Survey (SDSS; 
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Fig. 2. — (a) Mass surface density profiles for CGS Es with My < —21.3 mag. The individual profiles for the entire galaxy and for the combined inner+middle 
components are plotted in light grey and orange, respectively; their corresponding median profiles are plotted as darker solid lines, with the dashed lines denoting the 
associated 1 a uncertainty range of the median value (determined from bootstrap resampling), (b) Comparison of the median profile of the inner+middle components 
for CGS Es (orange solid line; dashed lines 1 a uncertainty) with the median profiles for massive galaxies at different redshift bins selected by their cumulative 
number density (iic = 1.4 X 10^ Mpc"'; Patel et al. 2013). The median profile (and associated uncertainty) for the 21 massive ETGsat 1.5 < z < 2.5 from Szomora 
et al. (2012) is also shown. Note that the radial scale is plotted as S' " in (a) but as linear in R in (b) to emphasize the outer regions. 



Guo et al. 2009, dashed line in Figure lb). As expected, high- 
Z ETGs, especially those with z > 1.5, fall systematically be- 
low the local relation. At fixed stellar mass, high-z systems 
are a factor ~ 4 - 5 more compact than their z = counterparts. 
While the inner component of the CGS galaxies is generally not 
massive enough to directly map onto high-z objects, the mid- 
dle component can be viewed as a slightly modified version of 
the compact quiescent galaxies. The mass-size relation for the 
middle component runs roughly parallel to, but is offset slightly 
below, that of local ETGs. It lies approximately in between the 
locus of high-z and low-z points. In dramatic contrast with the 
inner and middle components, the outer envelope displays sig- 
nificantly larger scatter on the Rg-M^, plane. While some of 
the scatter undoubtedly arises from measurement uncertainty, 
which, as indicated at the bottom-right comer of Figure 1(b), is 
considerable for extended, low-surface brightness features, we 
believe that most of the scatter is intrinsic. 

As discussed in Huang et al. (2013), the central core may 
have been reshaped more recently (Hopkins et al. 2009a, b, 
c) than the overall compact stellar body of interest, thus re- 
quiring a separate component to model well in order to sen- 
sitively extract the lower lying outer component. Regardless of 
the exact process that shaped the small central structure, it is a 
small perturbation to the overall structure of interest, hence we 
group the innermost two components (designated henceforth as 
inner+middle ) into a single entity. To measure the overall size, 
mass, and profile of the composite, we generated model images 
from best fit parameters that summed the innermost subcompo- 
nents, accounting for the PSF. Then, we fit a single-component 
model to extract the effective radius of this composite struc- 
ture. As the inner component has higher density than the middle 
one, we intuitively expect the effective radius of the combined 



structure to be smaller than that of the middle component alone. 
This is indeed the case, as Figure 1(b) shows. Surprisingly, the 
composite inner+middle structure (orange diamonds) not only 
forms a very tight mass-size relation with a slope similar to that 
of Guo et al. (2009) for SDSS ETGs, but it also Hes on top of 
the distribution of points for most high-z ETGs, especially those 
withz< 1.5. 

As advocated by van Dokkum et al. (2010), galaxies selected 
at constant cumulative number density provide a more physical 
connection between progenitors and descendants than those se- 
lected above a constant stellar mass. We therefore compare our 
observations with the sample of Patel et al. (2013), which is se- 
lected at a cumulative number density (n^ = 1.4 x lO"'^ Mpc""') 
that corresponds to massive ETGs (M* « 10^^-^ Mq) similar to 
the CGS sources at z ~ 0. Figure 1(b) shows that the median 
M* and R^ of our inner+middle components agree quite well 
with the properties of galaxies at 1.0 < z < 1.5 from Patel et 
al. (2013). Apart from global quantities such as mass and ef- 
fective radius, the surface density profiles themselves provide 
a more direct comparison of the structural connection between 
high-z ETGs and local Es. For this purpose, we convert the sur- 
face brightness profiles of the inner+middle components, after 
removing the point-spread function smearing, into stellar mass 
density profiles using the mass-to-light ratio adopted in Huang 
et al. (2013). Their median profile is shown in Figure 2(a), 
and in Figure 2(b) we compare it to the median profiles of the 
cumulative number density-selected galaxies in different red- 
shift bins from Patel et al. (2013). The median profile of the 
inner+middle components of local Es matches quite closely the 
profiles of massive galaxies at 1 .5 < z < 2.0 over most of their 
radial extent. Unlike the median M* and Rg, the median pro- 
file of 1.0 < z < 1.5 massive galaxies only marginally close to 
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the profile of inner+middle component. This seemingly differ- 
ence relates to the higher density of high-z objects at R < 1.0 
kpc. It is confirmed in the median density profile of the massive 
ETGs at 1 .5 < z < 2.5 from Szomoru et al. (2012), which show 
even higher central densities than the sample from Patel et al. 
at comparable redshifts. This apparent discrepancy may be due 
to different selection criteria employed in the two studies. 

The two-phase formation scenario of ETGs (Oser et al. 2010; 
Johansson et al. 2012) offers an appealing framework to inter- 
pret the above trends and to connect the distant and local pop- 
ulations of massive galaxies. Dissipative processes ("in situ" 
star formation, in the language of Oser et al. 2010) connected 
with early gas-rich events naturally lead to high central densi- 
ties. We suggest that the central, more compact structures of 
nearby Es — namely the "inner" or "middle" components, pos- 
sibly both — are the natural outcomes of this initial, dissipative 
phase. We propose that they are the direct remnants of the 
high-z massive quiescent galaxies. The combined luminosity 
(mass) fraction of the inner+middle components, / « 0.3-0.4 
are roughly comparable to the factor of ^ 2 increase necessary 
to grow the high-z objects, if they are the progenitors of lo- 
cal Es. Their overall physical scales are also consistent with 
their originally compact nature. In detail, the current evidence 
suggests that the high-z objects tend to have somewhat higher 
central densities than their local remnants. This is an impor- 
tant constraint that future theoretical models should try to re- 
produce. 

We equate the extended outer envelope with material ac- 
creted through dissipationless minor mergers ("ex situ" stars; 
Oser et al. 2010). The large intrinsic scatter on the R^-M^, 
plane is likely related to the stochastic nature by which the 
outer envelope is expected to have accumulated. The minor 
merger hypothesis has previously been questioned on the basis 
of the apparently small observed scatter of the mass-size rela- 
tion of nearby Es (Nipoti et al. 2012). For those two observa- 
tions to be consistent, dry mergers would have to preferentially 
add material to large radii, realized through many stochastic 
events spread over a protracted period. We propose that this 
history is represented by the extended, outer component identi- 
fied through our photometric decomposition. 

We end by briefly comparing our decomposition results with 
predictions from recent numerical simulations for the forma- 
tion of massive ETGs. After separating the final stellar system 
into its in situ (dissipative) and ex situ (accreted) constituents in 
their simulations, Oser et al. (2010) show that, in comparison 
to the in situ material, the typical size of the accreted compo- 
nent should correlate with the total stellar mass of the system. 
Furthermore, the stellar mass fraction of the accreted compo- 
nent may also correlate with the total stellar mass and effective 
radius of the whole system (Oser et al. 2010; Johansson et al. 
2012; Lackner et al. 2012). However, their numerical simu- 
lations produce a wide range of possible outcomes: the predic- 
tions for the variation of accreted fraction with total stellar mass 
in Oser et al. (2010) disagree strongly with those from the work 
of Lackner et al. (2012). 

We compare our analysis with the simulations to look for cor- 
relations in the accreted components. We assume that the outer 
and inner+middle components of our multi-component fits ap- 
proximate the "accreted" and "in situ" material in the theoret- 
ical models. Figure 3(a) shows that the effective radius of the 
outer component of nearby Es does indeed increase strongly 
with the total stellar mass of the system, much more so than 



either of the inner or middle components. While the phys- 
ical scales and the slope of the correlation are not the same 
as predicted in the simulations of Oser et al. , there is qual- 
itative agreement. However, the comparison between the ac- 
creted fraction, approximated as M^'^^VMlf ^', and effective ra- 
dius (Figures 3b) and total stellar mass (Figures 3c) do not 
show much, if any, correlation. Apart from the uncertainties 
of the predictions of current simulations, it is a difficult obser- 
vational challenge to isolate the accreted component via pho- 
tometric decomposition. The details of any decomposition are 
dependent on the choice of specific parameterized model com- 
ponents, which are not unique choices. Also, further galaxy 
merging would have diluted and partly erased the signatures. 

4. SUMMARY AND DISCUSSION 

Huang et al. (2013) show that the surface brightness distribu- 
tion of nearby Es can best be described by three photometrically 
distinct substructures, consisting of a compact {Re <^ 1 kpc) 
inner component, an intermediate-scale {Re « 2.5 kpc) middle 
component, and an extended {Re « 10 kpc) outer envelope. Mo- 
tivated by the recently proposed two-phase formation scenario 
for massive galaxies, here we argue that the inner and middle 
components identified in nearby Es jointly comprise a struc- 
ture that can be considered the evolved, local counterparts of 
the high-z compact massive galaxies. In this scenario, these 
compact structures derived from highly dissipative processes 
at high redshifts. Dissipationless minor mergers dominate the 
late-time evolution of ETGs by contributing to their dramatic 
size growth. We identify the outer, extended component of 
nearby Es with this late accretion phase. 

Several lines of investigation deserve further attention. In 
light of the photometric subcomponents uncovered in nearby 
Es, it is of interest to ask whether such substructure already 
exists in high-z ETGs, and, if so, whether their gradual build 
up can be detected as a function of time. Current studies treat 
high-z ETGs as single-component systems (e.g., Papovich et 
al. 2012; Szomoru et al. 2012), but the data may have suffi- 
cient quality to allow more complex analysis to reveal poten- 
tial substructure. Of particular interest is the mechanism by 
which the "core" structure in massive Es (e.g. Laueret al. 1995; 
Kormendy et al. 2009) arise. If scouring by binary supermas- 
sive black holes is necessary to make cores (e.g., Hopkins et 
al. 2009b, c), then major mergers are implicated. Similarly, if 
high-z massive galaxies contain disky or fast-rotating structures 
(van der Wei et al. 201 1), then they must be destroyed, presum- 
ably also through major mergers, if they are to evolve into local 
massive Es. How do major mergers fit into the framework of 
the two-phase formation scenario? 

Beyond imprints on the photometric structural properties, the 
two-phase formation scenario also implies differences in stellar 
population. Recent spectroscopic studies of the faint, outer re- 
gions of nearby Es already hint that they have lower metallicity 
and possibly older age than the inner parts (Coccato et al. 2010; 
Greene et al. 2012). Photometry is vastly simpler and less ex- 
pensive than spectroscopy, even if broad-band colors provide 
less robust constraints on stellar populations than spectra. A 
forthcoming paper in will make use of the multi-band CGS data 
to extract stellar population constraints from the color informa- 
tion of the subcomponents in nearby Es. 

We thank the anonymous referee for helpful comments and 
S. Patel for kindly providing us the data from his work. This 
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Fig. 3. — (a) Correlations between Re of different components and total stellar mass M^° . The symbols are the same as in Figure 1(a); the linear correlations 
from least-squares fitting are shown using the same color scheme as the points. The two right panels show the dependence between the mass fraction of the outer 
component and Re from single-component models (b) and MlJ" (c). The two lines in panel (c) represent the relation between accreted mass fraction (which we 
approximate with M°""='/M|,°'"') and M"''' from the simulations of Oser et al. (2010; long-dashed line) and Lackner et al. (2012; short dashed line). Note that for the 
purposes of this comparison we make use of the full sample 70 CGS Es with robust multi-component fits, not just the high-luminosity subsample. 
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